Dynamical properties of S = 1 bond-alternating Heisenberg chains in transverse 

magnetic fields 



Takahumi Suzuki and Sei-ichiro Suga 
Department of Applied Physics, Osaka University, Suita, Osaka 565-0871, Japan 

(Dated: February 2, 2008) 

We calculate dynamical structure factors of the S = 1 bond- alternating Heisenberg chain with a 
single-ion anisotropy in transverse magnetic fields, using a continued fraction method based on the 
Lanczos algorithm. In the Haldane-gap phase and the dimer phase, dynamical structure factors show 
characteristic field dependence. Possible interpretations are discussed in viewpoint of the field depen- 
dence of the excitation continuum. The numerical results are in qualitative agreement with recent 
results for inelastic neutron-scattering experiments on the S — 1 bond- alternating Heisenberg-chain 
compound Ni(C 9 D 2 4N 4 ) (N0 2 )C10 4 and the S = 1 Haldane-gap compound Ni(C 5 Di4N2)2N3(PF 6 ) 
in transverse magnetic fields. 
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I. INTRODUCTION 

5=1 bond-alternating Heisenberg chains have at- 
tracted a great amount of attention both theoretically 
and experimentally. The system shows a quantum phase 
transition between the Haldane-gap phase and the dimer 
phase depending on the bond-alternating ratio [|J . Ther- 
modynamic properties at the gapless point 0, LJ an d in 
the dimer phase 0, 0, IE were studied experimentally. 
The results were compared with theoretical results and 
good agreement was obtained 0, 0. The ele- 
mentary excitation was investigated by quantum Monte 
Carlo p and exact diagonalization jlj methods. In the 
Haldane-gap phase close to the uniform chain, the low- 
lying excitations are expected to be scattering states of 
the domain walls in the hidden antiferromagnetic order- 
ing 0. In the vicinity of the isolated dimer system, the 
low-lying excitations are well described by the S = 1 
magnon |9| . Dynamical structure factors in the Haldane- 
gap and dimer phases were calculated by the exact diago- 
nalization method 10]. In both phases, the distributions 
of the scattering intensity of the one-magnon mode and 
the excitation continuum were investigated. 

It is known that a feature of the Haldane-gap 
phase appears in transverse magnetic fields. In 
fact, the energy spectrum in transverse fields was 
measured for the anisotropic Haldane-gap compound 
Ni(C 2 D8N 8 )2(N0 2 )C104 (abbreviated to NENP) by in- 
elastic neutron-scattering experiments [ll| and the sat- 
isfying agreement between the theoretical and ex- 
perimental results was achieved. Many S — 1 bond- 
alternating Heisenberg-chains compounds possess notice- 
able single-ion anisotropy. Characteristic spin dynamics 
in the Haldane-gap and dimer phases may thus emerge 
in transverse magnetic fields. However, little knowledge 
had been obtained about dynamical properties of S = 1 
bond-alternating Heisenberg chains in transverse mag- 
netic fields. 

Quite recently, inelastic neutron-scattering 
experiments on the dimer-phase compound 
Ni(C 9 D 2 4N 4 )(N02)C104 (abbreviated to NTENP) 



in transverse magnetic fields (H) revealed fascinating 
aspects of dynamical properties In H = 0, the 

triplet magnon excitation is lifted into the lower S z = ±1 
branch and the higher S z = branch owing to a single- 
ion anisotropy. The observed scattering intensity of the 
higher S z = branch was anomalously weak at q = tt 
as compared with that in the anisotropic Haldane-gap 
compound Ni(C 5 Di4N 2 ) 2 N 3 (PF 6 ) (abbreviated to 
NDMAP), where the wave number is represented in the 
extended zone scheme. In H ^ 0, the degenerate lower 
branch is lifted. As H increases, the excitation energy 
of the highest (lowest) branch increased (decreased), 
while that of the middle branch was almost independent 
of H. In NTENP, the intensity of the highest branch 
decreased rapidly with increasing H and disappeared at 
H ~ 0.35i^ c , whereas the lower two branches survived 
up to H = H c . Note that H c is the critical field where 
the excitation gap closes. In H > H C1 there appeared 
only one gapped branch at q = tt which continued from 
the middle branch in H < H c . These findings are in 
contrast with those observed in NENP 11] and NDMAP 
|l4j under transverse magnetic fields: In H < H c three 
branches showed conspicuous intensity up to H = H c 
and in H > H c three gapped branches were observed. 
It is desirable to clarify the origin of the different spin 
dynamics between these systems. 

In this paper, we calculate the dynamical structure 
factor (DSF) of the S = 1 bond-alternating Heisenberg 
chain with a single-ion anisotropy in transverse mag- 
netic fields, using a continued- fraction method based 
on the Lanczos algorithm |l5|. By shifting the bond- 
alternating ratio systematically, dynamical properties in 
the Haldane-gap phase and the dimer phase are investi- 
gated. In Sec. II, we briefly summarize the method for 
the numerical calculation. In Sec. Ill, we show the re- 
sults for the DSF, turning our attention to the behavior 
at q — tt. In Sec. IV, the characteristic field dependence 
in the Haldane-gap phase and dimer phase is discussed in 
viewpoint of the field dependence of the excitation con- 
tinuum. The origin of different spin dynamics observed 
in NTENP and NDMAP is also discussed. Sec. V is 
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devoted to the summary. 



II. MODEL AND METHOD 



the gapless line than NDMAP 21]. In the dimer phase, 
NTENP is situated nearer the gapless line than NMAOP. 
In the next section, we show the results for S zz (q,uo) 
and S ± (q,cj) = S xx (q,uo) + S yy (q,u). 



Let us consider the 5 = 1 bond-alternating Heisenberg 
chain with a single-ion anisotropy in transverse magnetic 
fields described by the following Hamiltonian, 

H = J ^2 ($2i-i • S 2 i + aS 2 i • S 2 i+i) 

i 

+ dJ2(sD 2 -9»bhJ2s?, (i) 



where J > and a is a bond-alternating ratio. The 
periodic boundary condition is applied. The DFS can be 
expressed as |l5| 



S^(q,cu) 



1 



Im<W- 



1 



7T * Z - 

S^(q)C^(q,uj) 



(H = x,y,z), (2) 



where \^/ g ) is the eigenstate of TL with the lowest eigen- 
value E g , S% = (1/VN) J2j e iqj Sj with N being the to- 
tal number of spins, and z = Tiuo + irj + E g . The lattice 
constant between neighboring two spins is set to unity. 
Therefore, q = 0.57T corresponds to the boundary of the 
Brillouin zone in case of a ^ 1. We set gjiB = 1 and 
h = 1. The energy is measured in units of J. 

In the expression J2J), S^ v (q) is the static structure fac- 
tor and C^ LV {q 1 uo) is represented in the form of the con- 
tinued fraction, which can be calculated numerically by 
Lanczos algorithm. The total contribution of C fIU (q,uo) 
for a fixed q is normalized to unity, because the sum 
rule, S^ v {q) = J °° duoS^ u (q J w), has to be satisfied. In- 
stead of taking 77 — > +0, we set r] — 1.0 x 10 -2 . In 
finite-size systems, therefore, C flu (q,uo) consists of a fi- 
nite number of Lorentzians. The integrated value of each 
Lorentzian with respect to uo in given q is called as the 
residue [T^]. To discuss whether an excited state forms 
the isolated mode or the excitation continuum in the 
thermodynamic limit, the criterion proposed by Taka- 
hashi [la may be effective. According to this criterion, 
the residue of an excited state in the continuum tends 
to decrease with increasing the system size TV, while the 
residue for an isolated mode hardly depends on N. This 
method was successfully applied [l6| to explain spin dy- 
namics of NENP observed by inelastic neutron-scattering 
measurements fi% . 

Adopting the parameter sets (ol,D) = 
(1.0, 0.25), (0.8, 0.2), (0.45, 0.25), and (0.25,0.08), we 
calculate the DFS up to N = 20 spin systems in H = 
and TV = 16 spin systems in H 7^ 0. The parameter 
sets (a,D) = (1.0,0.25), (0.45,0.25), and (0.25,0.08) 
describe NDMAP [ll^NTENP Q, and the dimer-phase 
compound NMAOP pjj, respectively. According to the 
phase diagram [2(J, the parameter set (a,D) = (0.8,0.2) 
belongs to the Haldane-gap phase and is located nearer 



III. NUMERICAL RESULTS 

A. H = 

The DSF's in H = are shown in Fig. [I] Except 
for NDMAP, the results are shown in the extended zone 
scheme. The area of the circle is proportional to the scat- 
tering intensity. We turn our attention to the behavior 
at q = 7T. In (a, D) = (0.8, 0.2) and NTENP, the next- 
lowest-excited state in S zz (tt, uo) with relatively large in- 
tensity is located close to the lowest excited state, while 
in S- l (tt 1 uj) the next-lowest-excited state lies well above 
the lowest excited state. In NDMAP and NMAOP, there 
appear appreciable gaps between them in both S zz (tt,oj) 



and S- l (tt 
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FIG. 1: S zz (q,w) and S ± (q,u) at H = for N = 20. The 
area of the circle is proportional to the scattering intensity. 

In Fig. O we show the finite-size effects of the residues 
of the lowest excited states in S zz (q,uo) and S ± (q,uo). 
As shown in the insets, the residues of S zz (tt, uo) and 
S ± (tt j uo) in NDMAP and NMAOP show little size de- 
pendence, indicating that their lowest excited states form 
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FIG. 2: The size dependence of the residue in H = 0. Insets: 
The iV dependence of the residue at g = 7r. 



the isolated modes. In (a,D) = (0.8,0.2) and NTENP, 
on the contrary, the residues in S zz (tt, uj) decrease with 
increasing N, while the residues in S- l (tt,uj) scarcely de- 
pend on N. In spite of such appreciable size dependence, 
we have to be careful to decide the characteristic of the 
lowest excited state of S zz (tt, uj) in (a,D) = (0.8,0.2) 
and NTENP. Since their next-lowest-excited states are 
located close to the lowest excited states as shown in 
Fig. [I] the residues of the lowest excited states may be 
suffering from the higher energy states. 



The excitation energy and the intensity of the watched 
lower excited states are extrapolated to N — > oo. Their 
N dependence is shown in Fig. [3 We first evaluate 
the ratio of the intensity of the lowest excited state in 
S zz (tt j uj) to that in S ± (tt,uj) at H = 0. The lowest 
excited states in S zz (tt 1 uj) and S- l (tt 1 uj) at H = are 
the S z = and S z = ±1 branches, respectively. The 
resultant ratios are 0.35 for NDMAP, 0.29 for (a, D) = 
(0.8,0.2), 0.11 for NTENP, and 0.87 for NMAOP. Our 
results for NDMAP and NTENP are consistent with the 
experimental results [l3j . 
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FIG. 3: The intensity and the excitation energy (insets) of 
the lower excited states extrapolated to N — > oo in H = 0. 



B. H < H c 

The field dependence of the excitation energy and the 
intensity at q = ir is investigated in the same way. The 
results are summarized in Fig. 21 In magnetic fields, 
the lowest excited state in S ± (tt,uj) separates into two 
branches. As H increases, the lower branch shifts to 
the lower energy region and is softened at H C1 which re- 
sults in the quantum phase transition. The critical fields 
are evaluated as H c ~ 0.50 for NDMAP, H c ~ 0.40 for 
(a, D) = (0.8, 0.2), H c ~ 0.40 for NTENP, and H c ~ 0.65 
for NMAOP. Note that the gapped excitation energy is 
fitted well with while the softened mode is propor- 

tional to 1/N at H ~ H c 

We first discuss the behavior in S zz (ir, uj). As H in- 
creases, the lowest excited state in S zz (ir, uj) shifts to- 
wards the higher energy region, while the next-lowest- 
excited state shifts to the lower energy region. In 
NDMAP, the isolated mode never meets the excitation 
continuum even around H = H c . Therefore, the distinct 
intensity of the isolated mode appears in H < H c . In 
(a, D) = (0.8, 0.2), the isolated mode merges into the ex- 
citation continuum around H = 3.5 (< H c ), where the 
intensity of the isolated mode becomes so weak. Note 
that there are many excited states above the lower edge 
of the excitation continuum. However, their intensity is 
so small that they are invisible in the present scale. 

In NTENP and NMAOP, the intensity of the lowest 
excited states in S zz (tt, uj) diminishes with increasing H 
and disappears around H ~ 0.2 and ~ 0.4, respectively, 
where the lowest excited states cross the excitation con- 
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FIG. 4: The field dependence of the excitation energy and 
intensity for the isolated mode and the lower edge of the ex- 
citation continuum. In NTENP, the intensity of S zz (tt,lo) 
is enlarged by a factor of 5 as compared with the intensity 
in the other figures. The solid and broken lines denote the 
excitation energies for the isolated mode and lower edge of 
the excitation continuum, respectively. The critical fields are 
H c ~ 0.50 for NDMAP, H c ~ 0.40 for (a,D) = (0.8,0.2), 
H c ~ 0.40 for NTENP, and H c ~ 0.65 for NMAOP. 



tinua. After the crossing, the replaced lowest excited 
states reach uo ~ 0.4 at H ~ H c in NTENP. Judging 
from the results for (a, D) = (0.8, 0.2) and NTENP, at 
H = their lowest excited states in S zz (ir, uo) form the 
isolated modes. 

In 6' ± (7r, lj), on the other hand, the excitation continua 
of the four systems are located in the higher energy region 
up to H ~ H c . Accordingly, two isolated modes appear 
in H < H c in S ± (tv,uo). 



C. H> H c 

We next investigate dynamical properties in H > H c 
using the same method. In H > H c , the lowest-lying 
excitation becomes gapless, yielding the peak with the 
largest intensity in 5^(^,0) irrespective of H. We disre- 
gard the behavior of this peak. In Fig. we show the 



FIG. 5: The N dependence of the residue at q = tv in H > H c 
(left figures). H = 0.6 for NDMAP, H = 0.5 for (a,D) = 
(0.8, 0.2), H = 0.6 for NTENP, and H = 0.7 for NMAOP. The 
extrapolated intensity and the excitation energy in S zz (tv,cu) 
and S ± (ty,uj) are presented in the same figure (right figures). 
Insets: Extrapolation of the intensity (left figures) and the 
excitation energy (right figures) to N — > oo. 



results for NDMAP at H = 0.6, (a, D) = (0.8, 0.2) at 
H = 0.5, NTENP at H = 0.6, and NMAOP at H = 0.7. 
The residues of the peaks are shown in the left figures, 
and the extrapolated intensity and energy are presented 
in the right figures. The extrapolation of the intensity 
and the excitation energy are shown in the insets. 

In NDMAP, two appreciable peaks appear in S zz (ir, uo), 
while one peak appears in 6 ,± (7t,cj). Their residues are 
almost independent of N and the extrapolated intensity 
takes nonzero values. Therefore, these three peaks are 
from the isolated modes. Note that at H = H c the lower 
isolated mode in S zz (tt, uo) emerges at uo = 0, after the 
lower branch in S' ± (7r,cj) becomes softened. In (a,D) = 
(0.8,0.2), one appreciable peak appears in S zz (tt, uo) and 
6' ± (7r,a;), respectively. Their residues hardly depend on 
N and the extrapolated intensity takes nonzero values. 
Thus, these two peaks are from the isolated modes. Since 
the excitation continuum at H = 0.6 is located around 
oj = 0.8, only the lower branch of S zz (ir, uo) emerges. 

In NTENP, only one noticeable peak appears in 
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S- l (tt^uo). Its residue scarcely depends on TV, while the 
residue of the lowest excited state in S zz {tt 1 uo) decreases 
with increasing N. The results indicate that the low- 
est excited state in S ± (7t 1 cj) forms the isolated mode, 
while that in S zz (tt,uo) forms the lower edge of the ex- 
citation continuum. In fact, the intensity of the latter 
state extrapolated to N — > oo is much smaller than the 
extrapolated former one as shown in the inset. The lower 
edge of the excitation continuum of S zz (ir, uo) is evaluated 
as uo ~ 0.2 at H = 0.6, which makes the isolated mode 
in S zz (tt, uo) unstable in uo > 0.2. On the contrary, the 
excitation continuum of S- l (tt, uo) lies in uo > 1.2. Accord- 
ingly, in NTENP only one peak emerges in S- l (tt 1 uo). 

In NMAOP, one peak appears in S zz (tt, uo) and 
S- l (tt 1 uo) 1 respectively. Their residues hardly depend on 
TV, indicating that they are from the isolated modes. In 
H > H C1 the excitation continuum in S zz (ir, uo) shifts to 
the low energy region. At H — 0.7 its lower edge is lo- 
cated at uo rsj 0.2, which is close to the isolated mode. 
Therefore, as H increases in H > H c , the isolated mode 
in S zz (7r,uo) probably disappears. By contrast, the ex- 
citation continuum of S- l (tt,uo) lies in uo > 1.0, yielding 
the stable isolated mode in S- l (tt 1 uo) even in H > H c . 



IV. DISCUSSION 

S zz (q,uo) is active on the excitation process that con- 
serves the S z component of the ground state, while 
S ± (q^uo) is active on the excitation process that changes 
the S z component of the ground state by ±1. In H > H c , 
the ground state with (S x ) ^ and (S z ) = takes place 
instead of the singlet ground state in H < H c . In fact, we 
have confirmed that (S x ) =0-0.2 and (S z ) = O(10- 6 ) 
in H c < H < 0.7 for our four parameter sets. In H > H c , 
therefore, the excitations to the singlet state and to the 
triplet state with S z = component make main contribu- 
tions to the lower and higher isolated modes of S zz (tt, uo), 
respectively. 



A. Relation to inelastic neutron-scattering 
experiments 

We now discuss the observable DSF in the reduced 
zone scheme. In (a, D) = (0.8,0.2), the residues of 
S zz (q,uo) and S ± (q 1 uo) in q < OAtt show a typical pat- 
tern for the excitation continuum 16] as shown in Fig. 
12 The excitation energies of the isolated modes at q = tt 
are smaller than that around q ~ even in transverse 
fields. When < H < H C1 therefore, in the reduced 
zone scheme one isolated mode appears in S zz (0, uo) and 
two isolated modes appear in 5^(0,^). In H > H C1 one 
isolated mode emerges in S zz (0,uo) and 5^(0,^), respec- 
tively. 

In NTENP and NMAOP, the excitation energies take 
almost the same values at the symmetric wave numbers 
about q = 0.57T even in transverse fields. The same fea- 



ture was already known in the dimer phase for D = 
and H = [l(|. As shown in Fig. [21 the residues of 
S zz (q,uo) and S ± (q 1 uo) in q < OAtt show typical behavior 
for the isolated mode [l(|- Therefore, in the reduced zone 
scheme the lowest excited states of S zz (0, uo) and 5^(0, uo) 
form the isolated modes in H = 0. As H increases, in 
NTENP and NMAOP the intensity of the highest branch 
in S zz (0, uo) decreases and disappears at H ~ 0.2 and 
0.4, respectively. In H > H C1 only one isolated mode of 
5^(0, uo) emerges in NTENP, while the isolated modes of 
respective S zz ({),uo) and 5^(0, a;) emerge in NMAOP. 




FIG. 6: The field dependence of the excitation energies at 
q = TV. Full triangles denote S zz (tv,oo) and gray circles denote 
S ± (it,uj). On the solid and gray lines, appreciable scattering 
intensity appears. The broken lines for H c are guides to eyes. 

In inelastic neutron-scattering experiments on NTENP 
and NDMAP [H[3, the DSF's were observed as a su- 
perposition of S zz (q,uo) and S ± (q 1 uo). Using our numer- 
ical results, we summarize the corresponding results in 
Fig. |H1 On the solid and gray lines, the isolated modes 
possess appreciable intensity. In NTENP, the highest 
branch disappears at H ~ 0.5i^ c , which is larg er than 
the value ~ 0.35i^ c estimated experimentally [l3|. Above 
H C1 only one gapped branch emerges in NTENP, whereas 
three gapped branches appear in NDMAP. The field de- 
pendence of the DSF's shown in Fig. ^reproduce qualita- 
tivelythe experimental results for NTENP and NDMAP 

In NTENP 6] and NDMAP the phase transitions 
to the three-dimensional ordered states were observed by 
specific heat measurements in transverse fields H > H c 
at low temperatures. This phase transition is caused by 
the interchain interaction J'. In NTENP and NDMAP, 
their values are estimated as J' < 3 x 10 _3 J 24], and 
J' ~ 6 x 10 _4 J and J' < 1 x 10 _4 J depending on the 
direction 18], respectively. Judging from our numeri- 
cal results, the interchain interactions in NTENP and 
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NDMAP are so weak that they do not have serious ef- 
fects on dynamical properties in H > H c . 

B. Spin dynamics in the Haldane-gap phase and 
the dimer phase 

The characteristic dynamical properties are attributed 
to the different field dependence of the excitation contin- 
uum in S zz (tt^uj). As the system approaches the gapless 
line in both phases, the excitation continuum of S zz (tt, uj) 
easily shifts to the low energy region in weak fields. In 
the Haldane-gap phase this feature makes the highest iso- 
lated mode of (a, D) = (0.8, 0.2) unstable and invisible 
in H > H c , whereas in NDMAP three isolated modes ap- 
pear even in H > H c . When H > H Cl in the dimer phase, 
only one gapped mode of S^^^uj) appears in NTENP, 
while the lower isolated mode of S zz (tt, uj) emerges close 
to H = H c in addition to the isolated S- l (tt 1 uj) mode. 
As shown so far, the excitation continuum of S zz (ir, uj) 
in the dimer phase shifts to the lower energy region close 
to uj = as compared with that in the Haldane-gap 
phase. Such different field dependence of the excitation 
continuum may be intrinsic in the dimer phase and the 
Haldane-gap phase. 

V. SUMMARY 

We have investigated the DSF of the S = 1 bond- 
alternating Heisenberg chain with a single-ion anisotropy 



in transverse magnetic fields, using a continued fraction 
method based on the Lanczos algorithm. We have shown 
that the excitation continuum in S zz (q, uj) causes charac- 
teristic field dependence of the DSF in the Haldane-gap 
phase and the dimer phase. Our results well reproduce 
the different field dependence of experimental findings for 
NTENP and NDMAP. 
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